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Abstract 

The Generalized Pareto model is considered underlying 
model from which observables are to be predicted by 
Bayesian approach. The Type-II censored data from the 
model is considered for the Central Coverage Bayes 
prediction technique. Both the known and unknown cases of 
the parameters have been considered here. A simulation 
study also has been carried out for illustrating the 
performances of the procedures. 
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Introduction 

The Pareto model plays an important role in socio- 
economic studies. It is often used as a model for 
analysing areas including city population distribution, 
stock price fluctuation, oil field locations and military 
areas. It has been found to be suitable for 
approximating the right tails of distribution with 
positive skewness. The Pareto distribution and their 
close relatives provide a very flexible family of fat- 
tailed distributions, which may be used as a model for 
the income distribution of higher income group. The 
Pareto distribution has a decreasing failure rate, so it 
has often been used for model survival after some 
medical procedures (the ability to survive for a longer 
time appears to increase, the longer one survives after 
certain medical procedures). 

An important objective of a life-testing experiment is 
to predict the nature of future sample based on a 
current sample. Prediction of mean, the smallest or the 
largest observation in a future sample, has been a topic 
of interest and importance in the context of quality 
and reliability analysis. 

The objective of the present paper is to predict the 
nature of the future behaviour of the observation 
when sufficient information of the past and the 



present behaviour of an event or an observation is 
known or given for the generalized Pareto model. 
Under Type-II censoring scheme we obtain the 
Central Coverage Bayes prediction length of intervals. 
Both known and unknown case for the parameter has 
been considered. 

A good deal of literature is available on predictive 
inference for future failure distribution under different 
criterion. Few of those who have been extensively 
studied predictive inference for the future 
observations for the Pareto model are Arnold & Press 
(1989), Ouyang & Wu (1994), Mousa (2001), Soliman 
(2001), Nigm et al. (2003), Wu et al. (2004), Raqab et al. 
(2007). 

Model and Prior Distributions 

The probability density function and cumulative 
density function of the considered generalized Pareto 
model are given as 



f (x;o,e) = 



J_ 



F(x;o,6) = l- 



1- 



i 



;0<x<a,0>0, (2.1) 



(2.2) 



Here, is known as the shape parameter and c as 
the scale parameter. 

Suppose n items are put to test under model (2.1) 
without replacement and test terminates as soon as 

the r lh item fails (r < n) . If x p x 2 ,...,x r be the 

observed failure items for the first r components, then 
the likelihood function for first r failure 
items x(= x 1 ,x 2 ,...,x 1 .)is 



L(xie)oc6~ r o~ r exp 



I T 



(2.3) 
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where 



T,=-2>g 



a J 



+ (n-r) log 



1 — 



andT =log J 



f X ^ 



V 



When the shape parameter c is known, the scale 
parameter has a conjugate prior density and in 
present case it is taken as inverted Gamma 
distribution with the probability density function 



g(0) = 



p a e 



-p/e 



e a+1 r(a) 



; e>0,a>0,p >0. 



(2.4) 



We believe, as stated in Arnold and Press (1983) that 
from a Bayesian view point, there is clearly no way in 
which one can say that one prior is better than other. It 
is more frequently the case that, we select to restrict 
attention to a given flexible family of priors, and we 
choose one from that family, which seems to match 
best with our personal beliefs. Also the main reason 
for its general acceptability is the mathematical 
tractability resulting from the fact that the inverted 
gamma distribution is the conjugate prior for the scale 
parameter . Now, the corresponding posterior 
density for the parameter is obtained as 



*0) = 



(T r+ p) r+ "exp(-(T r+ p)/9) 
r(r + a) 9 



r + a + l 



(2.5) 



In most practical applications with proper informative 
information on the Pareto scale parameter, the prior 
variance is usually finite. Now we provide the prior 
information in the case when the shape parameter is 
unknown, which is more common in practice. It is 
known that in this case the Pareto distribution does 
not have a continuous conjugate joint prior 
distribution, although there exist a continuous-discrete 
joint prior distribution. The continuous component of 
this distribution is related to the scale parameter, and 
the discrete one is related to the shape parameter. 

In present case when the parameter a is unknown, 
the joint prior density for both the 
parameters and c is taken as 



gi(e,o) = 



P a e 



■p/e V i ^ 



e a+1 r(a) 



<pj 



6>0,0<o<^,(a,p, <p)>0. 



(2.6) 



The joint posterior density for the considered 
parameters and c is obtained as 



jr 1 (6,a) = 



0-(r + a + l) e -(T r + P)/9 e -T 



a r(r + a)c' 



e- T - 



a = 



(T r+ P) r+a 



da. 



(2.7) 



Central Coverage Bayes Prediction Interval 
(Parameter c is Known) 

Let x 1 ,x 2 ,...,x r be the first r observed failure items 
from a sample of size n under the Type-II censoring 
scheme from the model (2.1). If Y = ( y x , y 2 ,.„, y m ) 
be the second independent random sample of future 
observations from the same model. Then the Bayes 
predicative density of the future observation Y is 
denoted by h (y I x) and obtained by simplifying 

h(ylx) = Jf(y;o,e) jc(6) d6 



h(ylx) = 



(T r+ pr 



(3.1) 



(T r + p-log(G-y) + loga) r+a+1 ' 

In the context of Bayes prediction, we say that (/ P / 2 ) 
is a 100(1 — s)% prediction limits for the future 
random variable Y, if 



Pr(/j < Y</ 2 ) = l-e. 



(3.2) 



Here / , and / 2 are said to be lower and upper Bayes 
prediction limits for the random variable Y and 1 — S 
is called the confidence prediction coefficient. The 
Central coverage Bayes prediction lower and upper 
limits are obtain by solving following equality 



Pr(Y</ cl ) = ^ = Pr(Y>/ C2 ). 



(3.3) 



Here, I Cl and l C2 are said to be lower and upper 
Central coverage Bayes prediction limits for the 
random variable Y. Using (3.1) and (3.3), the lower 
and upper Bayes prediction limits of Y are obtained 

as 
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/ cl =o {l-e (T ' +p)(1 - El) } ; e,= 



a + eV l/(r + a) 



and 



Z t3 =o{l-ef r - + «( 1 -»)};8 2 = 



\ ^ J 



1-8 



(3.4) 



-l/(r + a) 



.(3.5) 



Therefore, the Central coverage Bayes prediction 
length of interval for the given future observation Y 
is 

^ c — ^C2 — ^ci • (3-6) 
Numerical Analysis 

To assess and study the properties of the Bayes 
prediction length of interval under the Central 
coverage the random samples are generated as follows: 

1. For the purposefully given values of the prior 
parameters a and [3, generate by using the 
prior density g ( ) . The considered values of 
prior parameters are (a,[3)= (02,0.50), 
(04,0.25) and (08,0.125). 

Here smaller values were consider for the prior 
parameters to reflect little prior information 
whereas larger reflect the higher prior information. 
Decisively, we kept the prior means the same as 
the original means, although we have seen (not 
reported here) that the overall behaviour remains 
the same. 

2. Using obtained in step (1), and the considered 
values of = 1.50,1.00,0.50 generates the 
10,000 random samples of size n(= 10) from the 
considered model (2.1). 

3. The selected set of censored sample size 

r = 04,06, 08,10 with the level of 
significance S = 99 %, 95 %, 90 % ; the Central 
coverage Bayes prediction length of intervals has 
been obtained and presented in Table 01. 

4. We observe from the table that the length of the 
interval tend to be closer as r increases or a 
decrease when other parametric values are fixed. 
Similar trend has also been seen when 
combination of the prior parameter increases. It is 
also noted that when confidence level decreases 
the length of intervals also decreases. 



Remark: 

In the case when the censored sample size r(= 10), 
the censoring criterion is reduces to the complete 
sample size criterion and hence the result are valid for 
complete sample case. 

Central Coverage Bayes Prediction Interval 
(Parameter a is Unknown) 

When both parameters are considered to be unknown, 
the Bayes predicative density for the future 
observation Y is denoted by h 1 ( y I X ) and 
obtained by simplifying 

h i( y lx ) = J \ f (y;a,0)7i 1 (0,a) d0 da 



a 9 



h 1 (ylx) = 



r + a 



(T r +p-log(a-y) + loga)- (r+a+1) 



a r e T ° (a-y) 



da. 



(5.1) 



Using (5.1) and (3.3), the lower and upper Bayes 
prediction limits for Y are does obtained by solving 
following equations 



V * J 



I a r (T r+ p) r + a 



1- 



1- (T^ 'lOg 



1 __ci 



>da 



and 
+ 



V * J 



J . cr r (T r+ p) r+a 



1- 



MTr+pO'log 



1- 



C2 



da 



For the future observation Y, the Central coverage 
Bayes prediction length of the interval is 



Ic — ^C2 ^ci 



Numerical Analysis 



(5.2) 



When both parameters considered as the random 
variable, a simulation study also has been carried out 
to study the properties of the Bayes prediction interval 
as follows: 
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1. For the similar set of the considered values of the 
prior parameters a and P , as taken in Section 4, 
the values of has been generated. 

2. For the given values of prior parameter 
(p (= 0.50, 1.00, 5.00), generate a random value 

of C from the uniform U ( 0, \l (p) . 

3. Using the above generated values of and G 
obtained in steps (1) & (2) we generates the 10,000 
random samples of size n(= 10) from the 
considered model (2.1). 

4. For the similar set of selected values for r and £ 
(Section 4), the Bayes prediction lengths of 
intervals have been obtain and presented them in 
the Table 02. 

5. It has been observed from the Table 02 that, all the 
properties of the Bayes prediction length of 
intervals are similar to the case when one 
parameter is known. 

Conclusions 

The Generalized Pareto model is considered here from 
the observables is to be predicted by Bayesian 
approach. The Type-II censored data is considered for 
the Central Coverage Bayes prediction technique 
under both known and unknown cases of the 
parameters. For illustrating the performances of the 
procedures a simulation study has been carried out 
under the purposefully choosing the values of the 
hyper parameters. 

The conjugate prior for the scale parameter is 
considered here on the basis of the general 
acceptability. In case when the shape parameter is 
unknown (a more common in practice), a continuous- 
discrete joint prior distribution is considered. 

For all the considered values it is observed that the 
central coverage Bayes prediction lengths of intervals 
are shorter. Also, the length tends to be shortest with 



decreasing in the sample size or the confidence level. 
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TABLE 1 BAYES PREDICTION LENGTH OF THE INTERVAL WHEN a IS KNOWN 





a = 1.50 


a = 0.50 


a = 0.25 


r 


a,[U 


99% 


95% 


90% 


99% 


95% 


90% 


99% 


95% 


90% 




2, 0.50 


5.2253 


3.6653 


1.8888 


4.5719 


1.0812 


0.5431 


0.8952 


0.5427 


0.2741 


04 


4, 0.25 


4.2984 


1.3686 


1.2375 


1.0558 


0.4996 


0.4139 


0.3688 


0.2891 


0.2580 




8, 0.125 


1.6773 


1.2530 


0.5750 


0.5042 


0.2613 


0.2314 


0.2639 


0.2465 


0.2140 




2, 0.50 


4.4177 


3.1275 


1.5639 


1.5260 


0.6585 


0.3658 


0.5901 


0.4772 


0.2438 


06 


4, 0.25 


2.1640 


1.3012 


0.6280 


0.5372 


0.4970 


0.3115 


0.2559 


0.1943 


0.1695 




8, 0.125 


0.9968 


0.7898 


0.4677 


0.4338 


0.1640 


0.1387 


0.1589 


0.1282 


0.1250 




2, 0.50 


3.2107 


1.2136 


1.0953 


0.8955 


0.4105 


0.3192 


0.3068 


0.2140 


0.1649 


08 


4, 0.25 


2.1192 


0.9216 


0.5288 


0.4574 


0.4039 


0.2396 


0.2253 


0.1657 


0.1574 




8, 0.125 


0.7699 


0.5345 


0.3577 


0.3570 


0.1570 


0.1334 


0.1563 


0.1260 


0.1114 




2, 0.50 


2.4252 


0.9805 


0.7441 


0.6816 


0.3855 


0.1954 


0.2397 


0.1887 


0.1618 


10 


4, 0.25 


1.9542 


0.8173 


0.5233 


0.3541 


0.2987 


0.1707 


0.1876 


0.1646 


0.1429 




8, 0.125 


0.7754 


0.3747 


0.3278 


0.2827 


0.1296 


0.1196 


0.1537 


0.1259 


0.1123 



TABLE 2 BAYES PREDICTION LENGTH OF THE INTERVAL WHEN a IS UNKNOWN 





cp = 0.50 


q> = 1.00 


cp = 5.00 


r 


a,[U 


99% 


95% 


90% 


99% 


95% 


90% 


99% 


95% 


90% 




2, 0.50 


4.2597 


3.7399 


1.9477 


3.6475 


1.1146 


0.5658 


0.8853 


0.5424 


0.2768 


04 


4, 0.25 


3.3237 


1.3965 


1.2761 


1.0771 


0.5151 


0.4312 


0.3647 


0.2889 


0.2606 




8, 0.125 


1.6937 


1.2785 


0.5929 


0.5143 


0.2694 


0.2411 


0.2609 


0.2463 


0.2162 




2, 0.50 


3.3908 


3.1517 


1.5925 


1.5378 


0.6704 


0.3763 


0.5764 


0.4710 


0.2431 


06 


4, 0.25 


2.1584 


1.3112 


0.6395 


0.5413 


0.5060 


0.3205 


0.2499 


0.1918 


0.1691 




8, 0.125 


0.9942 


0.7959 


0.4763 


0.4371 


0.1669 


0.1427 


0.1553 


0.1265 


0.1247 




2, 0.50 


3.2356 


1.2359 


1.1272 


0.9117 


0.4223 


0.3319 


0.3028 


0.2134 


0.1663 


08 


4, 0.25 


2.1357 


0.9385 


0.5442 


0.4657 


0.4155 


0.2492 


0.2224 


0.1652 


0.1587 




8, 0.125 


0.7758 


0.5444 


0.3681 


0.3635 


0.1615 


0.1387 


0.1543 


0.1257 


0.1123 




2, 0.50 


2.4142 


0.9862 


0.7563 


0.6855 


0.3917 


0.2007 


0.2337 


0.1858 


0.1610 


10 


4, 0.25 


1.9452 


0.8220 


0.5318 


0.3560 


0.3035 


0.1753 


0.1829 


0.1622 


0.1422 




8, 0.125 


0.7718 


0.3769 


0.3331 


0.2844 


0.1316 


0.1228 


0.1498 


0.1240 


0.1118 
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